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Phase behaviours of liquid crystals under external fields, conjugate to the nematic order and 
smectic order, are studied within the framework of mean field approximation developed by 
McMillan. It is found that phase diagrams, of temperature vs interaction parameter of smectic 
A order, show several topologically different types caused by the external fields. The influences 
of the field conjugate to the smectic A phase, which is fictitious field, are precisely discussed. 
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1. Introduction 

Electric and magnetic fields may influence on the liq- 
uid crystalline ordering and thus on the phase behaviour 
of the liquid crystal systems. These fields directly couple 
to the nematic order and influence on the nematic order- 
ing. Strong external fields may induce a nematic order 
even in the isotropic phase. 1,2 The external fields also 
slightly raise the isotropic-nematic (I-N) transition tem- 
perature. As the external field is increased, the nematic 
order in the nematic phase at the transition temperature 
decreases; together with the increase of the nematic order 
in the isotropic phase at the transition temperature, the 
discontinuity at the I-N transition is reduced. Eventually 
the discrete I-N transition vanishes under a sufficiently 
strong field, as shown in ref. 1 theoretically and in ref. 2 
experimentally. 

The influences of the external fields show a richer vari- 
ety and are much complex if we consider the smectic 
A order in addition to the nematic order. For exam- 
ple, a system which intrinsically shows direct isotropic- 
smectic A (I- A) transition under zero external field can 
exhibit the nematic phase between isotropic and smectic 
A phases under an external field. The existence of such a 
field-induced (non-spontaneous) nematic phase was pre- 
dicted theoretically in ref. 3 and found experimentally in 
ref. 4. The appearance of this non-spontaneous nematic 
phase is a consequence of the direct coupling between the 
field and the nematic order. The external fields also cou- 
ple to the smectic A order indirectly, through the direct 
coupling between nematic order and smectic A order. 
Thus the fields may influence the smectic order and may 
change the type of the nematic-smectic A (N-A) transi- 
tion. Actually, it is shown that a system which shows a 
first order N-A transition under zero external fields may 
undergo a second order N-A transition under an external 
field. 5 ' 6 

These studies have investigated the effects of the ex- 
ternal field conjugate to the nematic order, since the 
(squares of) magnetic and electric fields are the fields 
conjugate to the nematic order. 7 The external field con- 
jugate to the smectic order, which should induce a one- 
dimensional periodic density modulation in addition to 
the orientational order of the molecules, is not experi- 
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mentally available and not realistic field. (The electric 
field may play a role of the field conjugate to the smectic 
order by the effect of, so called, smectic electrostriction, 8 
which is the density change induced by the electric field. 
However, the effect of the electrostriction is negligibly 
weak except for the case where the coupling between ne- 
matic order and smectic order is weak. 8 ) 

The purpose of the present paper is to clarify the phase 
behaviour of the liquid crystals under the influence of not 
only the external field, h s , conjugate to the nematic order 
but also the external field, h a , conjugate to the smectic 
order. The field h a is fictitious and arbitrary field; how- 
ever, the application of such an external field is of interest 
from theoretical aspects since it completes the descrip- 
tion of the phase diagram. In addition, a field conjugate 
to an order parameter naturally arises as an effective 
field, as a result of the spatial inhomogeneity of the or- 
der parameter. Such an effective field is helpful to under- 
stand the phase transitions in inhomogeneous systems. 
In fact, one of the authors has disclosed the mechanism 
of phase transition of the freely suspended smectic films, 
by utilising the effective field conjugate to the smectic 
order parameter. 9 

In the present paper, we investigate the McMillan's 
liquid crystal model, 10 which has a parameter a show- 
ing a dimensionless interaction strength for the smectic 
A phase. By analysing order parameters and transition 
temperatures for a suitable range of a, we elucidate the 
phase behaviour of the system. The knowledge of the be- 
haviour of the homogeneous systems under the fictitious 
field is useful to investigate the inhomogeneous systems; 
the relation between the fictitious field and the effective 
field due to the spatial inhomogeneity is discussed in §4. 

2. McMillan Model 

We use a mean field approximation of liquid crystalline 
model introduced by McMillan. 10 Let s and a be order 
parameters of nematic order and smectic order, respec- 
tively. Then, the McMillan's one-body potential is ex- 
pressed as 



V(z,cose) = -v 



aa cos 



2?rz 



P 2 (cos6l), (1) 



where we assume the director is parallel to z-axis, and 
9 is the angle between the long axis of a molecule and 
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z-axis; P2{x) = 3(x 2 — l)/2 is the second order Legendre 
polynomial and d denotes the thickness of a single smec- 
tic layer. The second term on right hand side of eq. (1) 
denotes the smectic order interaction and parameter a 
is the smectic order interaction strength. The order pa- 
rameters are self-consistently determined from 



s = (-P2(cos0))/5 iS;(7 , 



a = { cos 



(2) 
(3) 



where the angular brackets denote the canonical ensem- 
ble average for one-particle: 



(A(6,z))^ 



1 



Z{p,pV s,(3V aa) 



I d9 sin9 / dzA(6,z) 
o Jo 



x exp <^ (3Vq 



s + acr cos 



2ttz 

~d~ 



P 2 (cos( 



(4) 



where Z denotes the one-particle partition function 



Z(f3,f3V s,pV a) 



= / d6»sin< 
Jo 

x exp < (3V 



V cxt = - 



dz 



[2tiz\ 
s + acr cos I — — I 



P 2 (cos6») 



(5) 



Under the external fields, the potential energy induced 
by the external fields is 

~~d 



h s P2(cos8) + h(j cos 



P 2 (cos 0) 



(6) 



where h s and h a are the external fields conjugate to, 
respectively, s and a. The mean field one-particle par- 
tition function is then Z(f3,f3h s + (iVos,[3h a + ftVoaa). 
For later convenience to obtain the free energy, we in- 
troduce symmetry breaking fields 12 r\ and (, which are 
conjugate to the nematic order and smectic order, re- 
spectively. The symmetry breaking fields should be zero 
at the thermal equilibrium. The one-particle partition 
function with symmetry breaking fields is Z([3,r) + [3h s + 
f3VoS, ( + f3h a + pVoaa). The self-consistent equations for 
order parameters are 



s = I(J3, i 1 + (ih s + f3V s, C + + pV aa), 
a = J(/3, r\ + {3h s + f3V s, ( + f3h a + f3V acr), 
where r], () and J(/3, 77, () are defined as 

J(P, V ,Q = ^-lnZ(p,r 1 ,0. 

Among the sets of solutions of cqs.(7), the thermo- 
dynamically stable set gives a minimum of the follow- 
ing function (the work done by the symmetry breaking 
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Fig. 1. Phase diagrams for (a)/i s = 0, h a = 0, (b)h s = 0.0106 > 
h( c \ h(j = 0, (c)h s = 0, h a = 0.003, and (d)h s = 0.0106, 
h a = 0.025. In (a), the isotropic, nematic and smectic A phases 
are indicated as I, N and A, respectively. The triple point and 
tricritical point are shown in (a) as X and +, respectively; these 
points for (a) are also shown on (b), (c) and (d) for reference. 
The circles indicate the tricritical point in (b), and the critical 
point in (c) and (d). 



fields): 

/3F((3, h s , h a ; s, a) 

= [ ds'r)(h s ,h (7 ;s',0)+ [ da'((h s ,h a ;s,cr'), (9) 
Jo Jo 

where rj(h s , h a \ s, a) and ((h S7 h a ;s,(T) arc the inverse 
solutions r\ and ( of the self-consistent equations (7). The 
minimum of F(f3, h s , h^; s, a) for given h s and h a is the 
thermodynamic free energy per molecule. We can rewrite 
eq.(9) by use of eqs.(7) as 

P (P{p, h a , K- s, a) - F((3, h s ,h a ; 0, 0)) 



f3V Q n/3V Q 2 
-—s + —aa 



-In 



Z{(3 7 ph s + /3V s, (3h a + pV aa) 



(10) 



where s and a are the solutions of eqs.(7) for given h s 
and h a with -q = ( = 0. 

3. Results 

Here, we show our numerical results obtained from the 
above formulation. For convenience of the following ob- 
servation, we first show the phase diagram of the sys- 
tems without external fields, which is originally shown 
in ref. 10. The a-T phase diagram is shown in Fig. 1(a). 
For a system with small a, i.e., a < 0.986, the system 
undergoes I-N transition at T = T m = 0.22019, and N- 
A transition at lower temperature Tna- The transition 
temperature Tna is a monotonically increasing function 
of a. For small a, the N-A transition is a second order 
transition, i.e., s and a are continuous but their differ- 
ential coefficients have discontinuity at the N-A transi- 
tion. If a and Tna are sufficiently large, i.e., a > a tc 
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and Tna > T tc , the N-A transition turns to first o: 
der transition. The tricritical values of a and Tna ai 
a tc = 0.703 and T tc = 0.1910. The existence of the tr 
critical point on N-A coexisting line was also pointed oi 
within the framework of the Landau-de Gennes theory 
The N-A coexisting line meets I-N coexisting line at tt 
triple point a tr i p i = 0.986, T tr i p i = 0.2202. A systei 
with a > a tr ipie exhibits direct I-A transition. 

We directly show that the independence of T IN on 
using the analogue of the Clapeyron relation for coe: 
isting lines of fluid phases. 11 As done in the Clapeyro 
relation we can write the slope of a coexisting line at 
point (a,Ti^(aj) as 

dT m (a) AF a 



da 



(1] 



where F is the thermodynamic free energy per molecu 
given by minimizing F(f3, h Sl h a \ s, a) with respect to 
and a. The factors on the right hand side of eq.(ll) are 
determined as AF a = <9F/<9a| ncm atic - dF/da\ isotm pic 
and AF T = dF/dT\ ncmat ic - <9T/<9T| iso tro P ic, where 
dF/da and dF/dT are evaluated at the point (a, Tin (a)) 
and at the phase indicated by the subscript. Then AF a — 
— A{Vq(t 2 /2) = since a is zero in both isotropic and 
nematic phases; AFt is the entropy change at the tran- 
sition and is nonzero. Therefore Tin is independent of 
a. 

In the following, we investigate the phase diagrams 
under the influence of external fields, h s and h a - We re- 
strict our investigations to ranges < h s < 0.03 and 
< h a < 0.06. We believe that the phase diagrams are 
qualitatively the same if the fields exceed the ranges, and 
that these upper bounds are sufficiently large to observe 
the influence of the fields. We also confirmed that, using 
our inhomogeneous model systems, the effective fields in- 
duced by the spatial inhomogeneity of the order param- 
eter do not exceed the ranges indicated above. The inho- 
mogeneous model systems we used are McMillan model 
systems sandwiched between strong homeotropic anchor- 
ing walls; we discuss these models and the effective fields 
in § 4. 

The phase diagrams of systems in nonzero external 
fields are shown in Figs.l (b), (c) and (d) for various h s 
and h a - These external fields change the phase diagram 
not only quantitatively but also qualitatively. 

We first discuss the effect of h s . If h s < h ( s c) = 0.0104, 
the application of h. s does not induce a topological change 
on the phase diagram, and increases both Tin and Tna 
linearly. We note that the triple point goes to high-a 
region by the application of h s . Thus, under the field 
h s , we have a parameter region in which the nematic 
phase appears, even if the systems do not exhibit the 
nematic phase under zero fields; i.e., the field induced 
non-spontaneous nematic phase 3 ' 4 appears also in the 
McMillan model. For h s > hi c) and h a = 0, the I-N co- 
existing line entirely vanishes and thus only N-A coexist- 
ing line remains (Fig. 1(b)). The fact that I-N transition 
vanishes under a strong external field h s was shown in 
rcf. 1 for a = 0, i.e., for Maier-Saupc model. 

The effect of h a is observed on the small a region even 
if h a is very small. The field h a removes the second or- 
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Fig. 2. The field h a dependence of critical values: (a)critical 
temperature AT NA (c) = T NA (c) Os,M - 7na (c) (h s , 0), 
and (b)critical parameter: AajMA.^- 1 = a NA^ c K^s , h a ) — 
aNA^C^sjO). The symbols o, □, o, and X correspond to h 3 = 0, 
0.011, 0.02, and 0.03, respectively. The curves for h s = and 
0.011 terminate at h a ~ 0.025 and ~ 0.04, respectively; this 
is because N-A critical points vanish for large h a since T^a^ 
increases and the critical points merge with the I-N transition 
line. 



der transition part on the N-A coexisting line (Fig. 1(c)). 
Thus the tricritical point turns to a critical point. As h a 
increases, N-A line moves to higher temperature region; 
I-N line is not sensitive to h a . If h a is strong enough, even 
if h s > h s ^ = 0.0104, a part of I-N line survives and it 
terminates at the I-N critical point (Fig. 1(d)). As h a 
increases, the temperature at I-N critical point Tin^ in- 
creases, and a value at I-N critical point ckin^ decreases, 
so that I-N line is elongated. The application of h s also 
increases both of Tn^ -* and aiN^, and thus shorten the 
I-N line. The aiN^ -* is sensitive to h s , so that the I-N 
line rapidly shrinks and vanishes for increasing h s . For 
h s > h s ^ and finite h„, there are paths where phases 
may change from isotropic to nematic, and nematic to 
smectic, without discontinuous transitions. Furthermore, 
I-N line rapidly shrinks as h s increases, and N-A line also 
rapidly shrinks as h a increases. Thus the bifurcated co- 
existing line becomes a single line which terminates at a 
critical point. 

We estimate the effect of h a on the N-A critical point 
more quantitatively. From our numerical calculations, 
the N-A critical temperature Tna^ is obtained to de- 
pend on h a as 

T NA ( c >(fo s ,/y)-JNA (c) (^0) _ ,2/5 no , 
T NA M(h„,0) ~ aK ' (i2) 

as shown in Fig. 2 (a). The ft. s -dependence of the coeffi- 
cient a is weak; the coefficients a are 0.423, 0.431, 0.458 
and 0.491 for h s = 0, 0.011, 0.020 and 0.030, respectively. 
The h s dependence of N-A critical point is approximately 
linear, i.e., for h a = 0, 



T NA (c) (frs,0)-?NA (c) (0,0) 
T N a< c H0,0) 



bh. 



(13) 



4 J. Phys. Soc. Jpn. 



Full Paper 



Author Name 



0.06 



0.04 



0.02 



1 1 1 

■ (a) 




A A 


I 


Ail 
i /, J 


i 





'.21 0.22 0.23 0.24 



0.21 0.22 0.23 0.24 0.25 
IcbT/Vo 



Fig. 3. The temperature-field (T-/i CT ) phase diagrams for (a)a = 

0. 975, h s = 0.01 < h{ c) , (b)a = 0.975, h s = 0.0106 > h{ c) , 
(c)a = 0.93, h s = 0.01, and (d)a = 0.93, h s = 0.0106. The labels 

1, N and A in (a) denote the isotropic, nematic and smectic A 
phases, respectively. 



where b = 1.88. Thus h s and h a dependence of Tna^ is 
completed as 

T NA ( c )(/i s ,/i ff )-r NA (c) (o,o) 



T NA (c) (0,0) 



bh s + o(l + &fc 8 )/# 5 . 

(14) 

If h s increases, the Tna^ increases up to I-N transition 
temperature, then N-A line vanishes. The /^-dependence 
of N-A critical value of a, ona*- ' , seems to be scalable for 
several h s as shown in Fig. 2(b). However, in the present 
stage, it is difficult to obtain a simple function form of 
scaled qna • 

In order to show clearly the effect of h a on the phase 
behaviour, we show h a -T phase diagrams in Fig. 3. In the 
system with a — 0.975, if the external field h s is weak, 
three phases are divided by first order transition lines 
(Fig. 3(a)). The transition temperature Tna increases 
almost linearly to h a . On the other hand, ^-dependence 
of Tin is weak. In fact, in the vicinity of h a = 0, Tin does 
not depend on h a . In analogy to the Clapeyron equation, 
we obtain 

dT m 



A{dF/dh a ) Aa 



where Aa and AS are the changes of, respectively, a 
and entropy per molecule at the I-N transition. Thus, 
dT m /dh a is zero since a is zero in both isotropic and 
nematic phases at h a = 0. By the application of h ai the 
smectic order is induced in isotropic and nematic phases, 
so that the change in a is finite and Tin becomes depen- 
dent on h a . If the field h s exceeds h s ^ c \ in the system 
with a = 0.975, the I-N transition turns to be contin- 
uous at h a ~ 0. But the application of h a recovers the 
I-N line thus the system again undergoes I-N transition 
(Fig. 3(b)). For h s < h s (c ^ and small a, e.g. a = 0.93, 
the N-A transition turns to be continuous under high h a \ 
thus N-A critical point appears (Fig. 3(c)). If h s exceeds 
h s ^ in a = 0.93 system, I-N transition at low h a be- 
comes a continuous change without transition. Then the 
system has paths from isotropic phase to smectic phase 



without discontinuous transition (Fig. 3(d)). 

4. Summary and Discussions 

We have studied the phase diagram of the McMillan 
model system under external fields h s , and h„. By the 
application of h a , second order transition line on N-A 
line at low a region vanishes and N-A transition turns to 
be smooth. The I-N line vanishes under an external field 
h s higher than h s ^ c \ however h a recovers the nematic 
order and I-N coexisting line. 

The fictitious field h a , which we mainly discussed in 
the present paper, is resulted from the spatial inhomo- 
geneity of the order parameter a. As an concrete ex- 
ample, let us consider a thin system sandwiched by two 
parallel walls and assume that the walls force liquid crys- 
talline molecules at wall surfaces to stick perpendicu- 
lar to the walls (homeotropic anchoring). In the follow- 
ing, we discretise the system into N layers parallel to 
the walls, and assume the nth layer has order parame- 
ters s(") and cr( n \ In the mean field theory, the inho- 
mogeneity due to walls changes the mean field contri- 
butions in the arguments of eqs.(7); e.g., Vqs turns to 
V#s (n) +Vi's {n -V +T/ " S ("+ 1 ) with constants V£ and y o " 
denoting intralayer and interlayer interactions, respec- 
tively. Then we may write the self-consistent equation 
as 

s («) = I(j3,0hW + f3V s {n \ I3h { a n) + PVqcutW), (16a) 

= J(p,0hW + pV s {n) , (3hi n) +pV aaW), (16b) 

where the and are defined by 

h{ n) = ^o"(s (n_1) - 2a<") + s ( ™ +1) ), (17a) 

= T/JV"- 1 ) - 2c» + a (n+ V), (17b) 

and the boundary conditions are s' 1 ^ = = 1 and 
1. In eqs.(16) and (17) we redefined 



r(l) 



Vq = Vq + 2Vq in order to emphasise the equivalence 
of eqs.(7) and (16) in the homogeneous limit. We sum- 
marised the effect of inhomogeneity of the order param- 
eters into h!f^ and h^. From the similarity between 
cqs. (7) and (16), we may call the variables and 
the effective fields corresponding to h s and h a , re- 
spectively. The stability condition for the inhomogeneous 
system is minimizing the work done by the symmetry 
breaking fields on the whole system, not the work on 
the each layer. Except for the difference in the stability 
conditions, each layer of the inhomogeneous thin system 
under effective fields may be viewed as a homogeneous 
bulk system under external fields which depend on po- 
sition and on temperature. From this point of view, we 
may analyse the thermal behaviour of the thin system by 
observing the locus of the fields on the bulk phase dia- 
gram of the field versus temperature plane. In fact, such 
an analysis is shown to be useful in order to uncover 
the mechanism of smectic A-C* phase transition, 9 and I- 
N phase transition of Maier-Saupe model. 13,14 For thin 
system version of the model we analysed in the present 
paper, we have investigated the phase transitions and 
thermal behaviour of physical quantities by utilising the 
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phase diagram obtained in the present paper, which will 
be reported in a near future. 15 
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